Anisotropic scattering of Bogoliubov excitations 
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We consider elementary excitations of an interacting Bose-Einstein condensate in the mean-field 
framework. As a building block for understanding the dynamics of systems comprising interaction 
and disorder, we study the scattering of Bogoliubov excitations by a single external impurity po- 
tential. A numerical integration of the Gross-Pitaevskii equation shows that the single-scattering 
amplitude has a marked angular anisotropy. By a saddle-point expansion of the hydrodynamic 
mean-field energy functional, we derive the relevant scattering amplitude including the crossover 
from sound-like to particle-like excitations. The very different scattering properties of these limiting 
cases are smoothly connected by an angular envelope function with a well-defined node of vanishing 
scattering amplitude. We find that the overall scattering is most efficient at the crossover from 
phonon-like to particle-like Bogoliubov excitations. 



Below a critical temperature Bose gases undergo a 
phase transition and form a Bose-Einstein condensate 
(BEC) [1-3]. It is a long-standing question how such 
a condensate is influenced by the competition between 
inter-atomic interactions on the one hand and external 
disorder on the other [4, 5]. Generically, both the ground- 
state phase diagram and non-equilibrium features depend 
crucially on the presence of soft modes and their proper- 
ties [6, 7]. At low temperatures, the relevant excitations 
of a BEC are Bogoliubov excitations [8] with collective 
properties due to the repulsive inter-atomic interaction. 
Their dispersion relation interpolates between the col- 
lective sound-wave and the single-particle regimes. In 
gaseous BEC, these Bogoliubov excitations can be exper- 
imentally created and analysed using Bragg spectroscopy 
[9-12]. ' 

In this article we focus on the controlled 2D scattering 
of Bogoliubov excitations by a single elementary impu- 
rity. We present numerical evidence for highly anisotro- 
pic scattering that features a characteristic node. By a 
variational treatment of the quantum hydrodynamic en- 
ergy functional, we derive the effective Hamiltonian for 
impurity scattering and obtain analytical expressions for 
the scattering amplitude and the position of the node. 
We expect our results to be useful in the future for under- 
standing the quantum-transport dynamics in disordered 
media, which builds on repeated single-scattering events 
[13]. 



SIMULATION 

We start our investigation by numerically simulating 
the scattering of Bogoliubov excitations by an impurity 
potential in a two-dimensional BEC, see Fig. 1(a), un- 
der periodic boundary conditions that mimic the very 
shallow trap required for an experimental realisation. At 
temperatures much lower than the critical temperature 
[1-3], a mean-field description in terms of the macro- 
scopically occupied wave function ^(r,t) is appropri- 
ate. First, we obtain the static ground-state density 
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Figure 1: (a) A blue-detuned laser beam focused perpen- 
dicularly through a 2D BEC provides an impurity potential 
V(r) = Vexp(— r 2 /ro) that scatters an incident Bogoliubov 
wave with wave vector k into modes k! at scattering angle 6. 
(b) Density plot \5 1 i/ k i j 2 of the scattered state obtained by nu- 
merical integration of the GP equation. The components of 
the scattered wave are distributed on the elastic circle |fe'| — k 
with characteristic nodes at ±#o- For the numerical param- 
eters V = 0.25^t, fc£ = 1, kro = 0.5, scattering is suppressed 
around 9o m tt/3. 



n o(r) = l^o | 2 °f the condensate in presence of an 
impurity potential V(r) = Vexp(— r 2 /rg) by imaginary 
time propagation of the Gross-Pitaevskii (GP) equation 
[1, 14]. Then a plane wave Bogoliubov excitation with 
wave vector k is superimposed on this ground state so- 
lution by a suitable choice of initial conditions at time 
t — 0. Then the ensuing time evolution according to the 
full GP equation is calculated. During the simulation the 
Bogoliubov wave moves forward and is scattered at the 
impurity. 

In order to analyse the scattered state, we Fourier 
transform a snapshot of the deviation S^(r, t) = ^(r, t) — 
^o(r) from the ground state around the impurity. 
Fig. 1(b) shows the resulting momentum density IJ'I'fc'l 2 . 
The scattering is essentially elastic, with the components 
of the scattered wave distributed on the circle \k \ = k. 
Surprisingly at first sight, scattering is suppressed at two 
symmetric angles ±8q with respect to the forward direc- 
tion. In the deep sound-wave regime fc£ <C 1 (where £ 



2 



is the healing length), we find these nodes at 9q — ► 7r/2 
resulting in a dipole scattering (p-wave) characteristic. 
For large values of fc£, when the Bogoliubov excitations 
are particle-like, the nodes shift to the forward direction, 
#o — > 0. This intriguing anisotropic scattering of Bogo- 
liubov excitations is a signature of the intricate crossover 
from single-particle to collective excitations in interacting 
Bose-Einstein condensates. 



LIMITING CASES 

Before tackling the full quantum hydrodynamical 
problem, we find it instructive to discuss the two lim- 
iting cases of pure particle-like excitations 3> 1) and 
pure sound- like excitations <C 1), where the expected 
angular distribution can be derived from elementary con- 
siderations. 

In the single-particle part of the Bogoliubov spectrum, 
excitations arc plane matter waves with dispersion rela- 
tion e° = H 2 k 2 /2m. The amplitude of a single-scattering 
process k i— > k' = k + q is proportional to the Fourier 
component V q = L~ d / 2 j d d rV{r)e~ lqr of the scattering 
potential. If the potential V(r) varies on a characteristic 
length ro , the scattering may be anisotropic if the wave 
can resolve this structure, qr n > 1 [15, 16]. In the op- 
posite case of a small obstacle such that qr -c 1, also 
known as the s-wave scattering regime, the scattering 
amplitude is simply proportional to Vq and can therefore 
only be isotropic. 

Quite on the contrary, we expect a very anisotropic 
scattering amplitude, proportional to Vq cos 9, in the deep 
sound-wave part of the Bogoliubov spectrum. Indeed, in 
the Thomas-Fermi regime, where the healing length £ is 
much smaller than the scale ro of potential variations, 
the condensate ground-state density follows the poten- 
tial: n^p(r) = [/j, — V(r)] j g. (Here as in the following, 
we assume that the impurity potential is always smaller 
than the chemical potential /i, such that n > every- 
where.) Excitations of the superfluid ground state in 
the regime fc( « 1 are longitudinal sound waves with 
density fluctuations Sn = n — «tf and phase fluctua- 
tions Sip. Importantly, the phase is the potential for 
the local superfluid velocity v s = (h/m)'VS(p. The su- 
perfluid hydrodynamics is determined by the continuity 
equation 9 t n + V • (nv s ) = on the one hand, and 
by the Euler equation for an ideal compressible fluid, 
m [dtV s + (v s ■ V)d s ] = — V(gn + V), on the other. To 
linear order in Sn and Sip, these two equations can be 
combined to a single wave equation 



m [c 2 V 2 - a 2 ] Sn = V ■ (V(r)VSr 



(1) 



where the sound velocity appears as c = \J [x/m. The 
excitations of a homogeneous fluid (V(r) = 0) are 
plane sound waves with linear dispersion u>k = ck. 



The gradient-potential operator on the right-hand side 
then causes scattering with an amplitude proportional to 
— (fc • k')Vk'-k- Hence, the potential component Vy-h, 
which must appear in all cases to satisfy momentum con- 
servation, is multiplied with a dipole (or p-wave) charac- 
teristic A{9) = — cos 9. This scattering cross-section with 
a node at 9q = ±ir/2 can be understood, in the frame of 
reference where the local fluid velocity is zero, as the 
dipole radiation pattern of an impurity that oscillates to 
and fro, quite similar to the case of classical sound waves 
scattered by an impenetrable obstacle [17]. 

By continuity, there should be a smooth crossover 
from the sound-wave behaviour A(9) = —cos 9 to the 
single-particle case A{&) — 1 as the excitation wave vec- 
tor k explores the Bogoliubov dispersion relation = 
Vefc(eg + 2/x). In the following, we will derive the corre- 
sponding analytical expressions for the relevant envelope 
function A(k£, 9) and the position of the scattering node, 

0o 



VARIATIONAL THEORY 

Since the scattering node is clearly present in the hy- 
drodynamic regime, we choose the density-phase repre- 
sentation ^(r, t) = n(r, t) exp{i(p(r, £)} and start with 
the grand canonical energy functional [18] 
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The chemical potential fi determines the total number of 
particles and introduces the healing length £ = h/y/2mfi 
as the length scale on which the condensate can respond 
to a spatial perturbation. The interaction constant g > 
stabilises the superfluid behaviour of the condensate. 
The external potential V(r) shall describe the local im- 
purity, with the influence of the very shallow trap in the 
centre of the BEC being negligible. In order to describe 
the dynamics of Bogoliubov excitations in the presence 
of an impurity potential, we use a four-step procedure (i- 
iv), equivalent in spirit to [18], but with results somewhat 
more useful in the present context. 

(i) The condensate ground state density n (r) and 
phase ipo are determined in the presence of the exter- 
nal potential V(r) as the saddle-point solution of the 
mean- field energy functional (2), 
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One finds that the kinetic energy is always minimised by 
a spatially homogeneous phase V^o = 0, i.e., absence 
of superfluid flow. The ground-state density no(r) as 
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function of V{r) solves the stationary equation 

h 2 (v 2 v^) 



2m 



+ gn = n - V. 



(4) 



(ii) Density fluctuations 6n(r,t) = n(r,t) — no(r) and 
phase fluctuations 6(p(r,t) = <p(r,t) — ipo are conjugate 
variables that obey the coupled equations of motion 
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(5) 



The relevant energy functional F is obtained by a 
quadratic expansion E = E + F[Sn, Sip] around the sad- 
dle point and reads 



(v 2 v^) 



3/2 



Sn 2 



+ gSn 2 }. (6) 
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In this formulation, the external impurity potential V(r) 
affects the fluctuations only through its imprint (via (4)) 
on the ground-state density n (r), and this visibly in a 
highly nonlinear manner. 

(iii) Since we wish to calculate the scattering amplitude 
to linear order in V, it also suffices to know the ground- 
state density to the same order. As shown in [19], by 
linearising (4) for small deviations from the bulk density 
n>oc = fJ'/g-, the condensate density no{r) can be written 
in Thomas-Fermi form 



n Q {r) 



[1 



;(r)]. 



(7) 



Here, the smoothed dimensionless potential v(r) = 
V{r)/n is a convolution of the bare potential V(r) by 
a Green's function with a very simple form in fc-spacc: 



V k /fi 



l + k 2 e/2 



(8) 



This formula, derived in different notations already some 
time ago (eq. (11) in [18]), shows that Fourier compo- 
nents of the effective potential with k ^S> l/£ are sup- 
pressed. In other words, the condensate does not follow 
features of the potential varying on a length scale shorter 
than the healing length £. 

Using the smooth ground-state density (7) in the en- 
ergy functional (6) and developing all terms to linear or- 
der in v(r), we can write F = F^ + F^ as the sum of 
two terms: the energy of excitations of the homogeneous 
bulk condensate with density rioc = n/g, 



F 



(0) 



d d r 



2m 



4n„o 



9 -Sn 2 
2 
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and a perturbation where the smoothed impurity poten- 
tial couples linearly to several gradient terms: 
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(10) 

This scattering term F^> is quadratic in the fluctuations 
and linear in the external potential and thus goes beyond 
Huang and Meng's theory [20]. 

(iv) The free-space contribution F' ' can be diago- 
naliscd by going into Fourier modes followed by the Bo- 
goliubov transformation 



iSifk 
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(11) 



Choosing a& = efe/ (2e°.rc.oo) in terms of the single-particle 
energy = h 2 k 2 /2m and the Bogoliubov dispersion = 
V e k( e k + then indeed gives 
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(12) 



At this point, these Bogoliubov excitations could be con- 
veniently quantised by imposing canonical commutation 
relations, which is not needed for the present purpose 
such that we continue to treat the 7^ as complex field 
amplitudes. 

Upon Fourier-Bogoliubov transforming, the impurity 
scattering contribution (10) acquires the structure 
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7,7 



L d/2 



7fc , 



(13) 
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plus terms containing products ^k'lk and 7^7^ which 
can be disregarded for scattering to linear order in the 
impurity potential V. We find that the elastic scat- 
tering amplitude as function of the on-shell momenta 
|fe| = |fe I = k and the scattering angle 9 = Z(fc,fe') 
writes 

W(k,9) :=Wk'k\ k , =k = ^A(kt9)V(k,9). (14) 

Here, the expected potential factor V(k, 6) — V q at q = 
\k' — k\ = 2fcsin(£>/2) is completely factorised from a 
remarkably simple angular envelope, 



A(kC,( 
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(15) 



This angular envelope, drawn as a polar plot for sev- 
eral values of fc£ in Fig. 2(a), describes the smooth tran- 
sition from sound-wave to free-particle scattering as a 
function of reduced momentum fc£. In the deep sound- 
wave regime fc£ — > 0, A = — cos(#) reproduces the dipole 
radiation pattern predicted by the hydrodynamic equa- 
tion (1). There is always a sign change between forward 
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Figure 2: (a) Polar plot {A(k£,8)] 2 of the angular envelope 
function (15) for fc£ = 0.2 (red), 0.5, 1, 2, 5 (violet). The enve- 
lope is close to a dipole radiation (p-wave) pattern for sound 
waves fc£ <C 1, and tends to an isotropic (s-wave) pattern 
for single-particle excitations fc£ 3> 1. In the intermedi- 
ate regime, backscattering is favoured over forward scatter- 
ing, (b) Density plot of Im[7 fe /], the imaginary part of the 
Bogoliubov-transformed amplitude = uy&^y +v k /8't>*_ k , 
from the numerical simulation of the GP equation. This fig- 
ure is more distinct than Fig. 1(b), because the interference of 
Bogoliubov waves with ±fc' is eliminated, and clearly shows 
the sign change across the scattering node 9q. Numerical pa- 
rameters as in Fig. 1(b). 




Figure 3: Coarse-grained elastic scattering amplitude (17) in 
units of / = 7oVb/(8£/i) for different values of at fixed 
potential radius kro — 0.5. Symbols: Results from numerical 
integration of the full GP equation. Solid curves: analytical 
prediction (17). With increasing k, the node moves to the left, 
according to (16). The overall amplitude has a maximum at 
fc£ w 1. 



and backward scattering since A(kt;,ir) = 1 = — A(k£,0) 
holds independently of fc£. From (15), the resulting node 
of vanishing scattering amplitude is found to be at 

k 2 e , , 



In the particle regime fc£ 3> 1, the nodes ±6o shift to 
the forward direction such that A(k£, 9) converges point- 
wisely to the isotropic single-particle envelope A = 1 
as fc£ — > co. Finally, when the healing length £ be- 
comes larger than the system size L, the node angle 
9q ps s/2/k£ becomes smaller than the angular fc-space 
resolution 1/kL. Then, the last contribution with nega- 
tive A(fc£, 9) is the forward scattering element, which can 
be absorbed by shifting the origin of the single-particle 
energy e k , and we recover the Hamiltonian for the poten- 
tial scattering of free matter waves [15, 16]. 

BACK TO THE SIMULATION 

In order to confront these predictions with the nu- 
merical results, we first of all Bogoliubov-transform the 
numerically calculated momentum amplitude, j k > = 
Uk'^k'+^k'S^*^ , with the usual coefficients u k = (efe + 

el)/{2^ k ) and v k = (e fc - e°)/(2V^i£). In the cor- 
responding plot of Imfy^'], the imaginary part of the 
Bogoliubov-transformed amplitude, shown in Fig. 2(b), 
one can clearly see the amplitude sign change across the 
scattering node. 

For a quantitative comparison with the numerical 
simulation, let now a single Bogoliubov excitation 



7k (t) = 70 exp{— iekt/h} with wave vector k be scat- 
tered by the impurity potential as described by the 
effective Hamiltonian (12) and (13). In linear re- 
sponse (equivalent to the Born approximation of the 
corresponding quantum problem), the scattered ampli- 
tude %?{t) = 7^,7 exp{-ie k t/H} is given by 7^ = 
7 G (e fe , k')L- d l 2 W k * k , where G (e, k) = [e-e k + iOf 1 
designates the retarded Green function of the free propa- 
gation described by (12). For a quantitative comparison 
with the numerical simulation in the finite system, one 
can coarse-grain over a fc-space area fcA#Afc around a 
point on the elastic circle at a given angle 9. The ana- 
lytical prediction for the coarse-grained imaginary part 
then reads 

lm[ 1 ^(9)} = -n l0 ^^W(k,9) (17) 

in terms of the density of states per unit area, Po( e ) = 
— 1/(ttL 2 ) J2 k Im[Go(e, k)]. Coarse-graining similarly the 
numerical results on the elastic circle, we can plot to- 
gether both amplitudes as function of the scattering angle 
9 at various wave vectors fc£, see Fig. 3. The agreement is 
clearly very good, with residual numerical scatter around 
the analytical curves due to transients and boundary ef- 
fects. 

The overall magnitude of the scattering amplitude in 
Fig. 3 first grows and then decreases as fc£ crosses from 
the sound-wave to the particle-regime, with most efficient 
scattering for wave vectors fc£ ks 1. This behaviour re- 
sults from two competing scalings: the Bogoliubov scat- 
tering amplitude W(k, 9) cx e°/efc is proportional to k 
for k( <C 1 and saturates to a constant for fc£ ^> 1. 
The factor po{c k )/k on the other hand is proportional 
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to the inverse group velocity (dk/dek) that behaves like 
the constant c _1 for sound waves k( < 1 and decreases 
as fc _1 for particles k£ ^> 1. The product of both con- 
tributions in (17) therefore has limiting behaviour k and 
k~ 1 , respectively, with a scattering maximum around the 
crossover fc£ w 1 from phonons to particles. 

Note that Fig. 2(b) is much clearer than Fig. 1(b), 
because there Bogoliubov excitations with opposite ±fc' 
interfere in the wave function densities \8^>y | 2 . One may 
wonder why the superposition of nodes stemming from 
opposite wave vectors still gives a density dip as clear as 
in Fig. 1(b). In fact, in the single-particle case fc£ ^ 1 the 
ratio Vk /uk tends to zero such that only the node of one 
component is observed, whereas for sound waves fc£ <IC 1 
both components contribute equally, but now with sym- 
metric nodes at ±-| that superpose exactly. This node ro- 
bustness should facilitate the experimental observation. 



single-scattering amplitude factorises into the impurity 
part and an angular envelope that describes the continu- 
ous transition from wave to particle behaviour as function 
of excitation momentum. Due to this factorisation, the 
theory is independent of the actual shape of the potential 
and can also be employed to describe disordered systems. 
In particular, we plan to generalise recent results on the 
disorder-induced localisation of Bogoliubov excitations in 
ID [23, 24] to higher dimensions. Especially the 2D case 
promises to be interesting, because this is the lower crit- 
ical dimension for the Anderson model of noninteracting 
particles in a random potential [25]. 

Financial support from DFG, BFHZ-CCUFB, and 
DAAD is gratefully acknowledged. We thank V. Gurarie 
for explaining the dipole scattering pattern with the hy- 
drodynamic formulation that proved to be very fruitful. 



EXPERIMENTAL REALIS ABILITY 

We propose our theoretical predictions to be experi- 
mentally tested following the numerical setup: A mod- 
erately strong, blue-detuned laser is focused perpendic- 
ularly through a 2D condensate [2, 21] without deplet- 
ing the condensate entirely. Then Bogoliubov excitations 
arc imprinted optically [9-12] and observed in a sub- 
sequent time-of- flight measurement at time to- If nec- 
essary, the sensitivity to certain fe-components can be 
greatly improved using Bragg spectroscopy [10, 22]. In 
both cases, the total momentum distribution l^fc^o)] 2 = 
l^fc + <5\I/ fc (io ) 1 2 a t time to is accessible. This is an os- 
cillating quantity since the inverse Bogoliubov transfor- 
mation S^kit) = Ufc7fc(t) — Wfe7lfc(0 superposes compo- 
nents of Bogoliubov waves scattered into opposite direc- 
tions with conjugate phases. The excitations S^k(t) live 
on the stationary background of the impurity-deformed 
condensate ground state ^ oc Vk given by (7) such that 
the time-of-fhght density would read |* fe (t )| 2 = (*fc) 2 + 
2^ k Rc [S 1 ^ k{to)] + \^ k\ 2 ■ If a time average is performed, 
the linear oscillating term drops out, and the density 
l^fcl 2 can be extracted by subtracting the ground-state 
density. Choosing an appropriate measurement time to 
can reveal the amplitude Re [<5 v I / fc(io)] + 0(8^ 2 ) on the 
smooth background ^ with a better signal-to-noise ra- 
tio. The experiment could notably test the limits of va- 
lidity of the weak-scattering linearisation and more gen- 
erally the breakdown of mean-field behaviour. 

CONCLUSIONS 

A variational treatment has allowed us to derive a sim- 
ple impurity-scattering Hamiltonian that governs the dy- 
namics of Bogoliubov excitations in Bose supcrfluids in 
presence of a weak external potential. Remarkably, the 
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